Electronic State and Magnetic Susceptibility in Orbitally Degenerate (J = 5/2) 
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Magnetic susceptibility in a heavy fermion system is com- 
posed of the Pauli term (xp) an d the Van-Vleck term (xv)- 
The latter comes from the interband excitation, where /- 
orbital degeneracy is essential. In this work, we study xp and 
Xv in the orbitally degenerate (J = 5/2) periodic Anderson 
model for both the metallic and insulating cases. The effect of 
the correlation between /-electrons is investigated using the 
self-consistent second-order perturbation theory. The main 
results are as follows, (i) Sixfold degenerate model: both xp 
and xv are enhanced by a factor of 1/z (z is renormalization 
constant), (ii) Nondegenerate model: only xp is enhanced 
by 1/z. Thus, orbital degeneracy is indispensable for the 
enhancement of xv- Moreover, orbital degeneracy reduces 
the Wilson ratio and stabilizes the nonmagnetic Fermi liquid 
state. 
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In general, the Ce-compound heavy fermion (HF) sys- 
tem possesses one electron per Ce-atom. Due to strong 
l-s coupling, the lowest /^eigenstate has the angular mo- 
mentum J = 5/2. In real materials, the sixfold states in 
J = 5/2 are split into three Kramers doublets (KDs) by 
a weak crystalline electric field (CEF). The higher KDs 
strongly affect the electronic structure because /-level 
splitting is, at most, ~ O(10 2 )K. In HF systems, orbital 
degeneracy (or orbital freedom) has an essential effect 
on many physical properties. For example, it causes the 
anomalous Hall effect, or a quadrapole long range order. 
Another interesting feature of the orbitally degenerate 
systems is the existence of the Van-Vleck susceptibility 
Xv, which is brought about by magnetic interband exci- 
tation. 0-^] xv should bring about large residual mag- 
netic susceptibility in many Kondo insulators (CeNiSn, 
Ce3Bi4Pt3, etc), || or in superconducting compounds 
below T c (UPt3, UBei3, etc). Thus, orbital freedom 
causes various nontrivial phenomena in HF compounds, 
and should be studied in earnest. In addition, orbital 
freedom also exists in many d-electron systems (V- or 
Mn-oxide compounds, etc). 

In this work, we study the orbitally degenerate (J = 
5/2) periodic Anderson model (J = 5/2 PAM), which 
is the appropriate model for Ce-compounds. In this or- 
bitally degenerate model, the magnetic susceptibility x 
is composed of the Pauli term (xp) and the Van-Vleck 



term (xv)- It is proved that xv appears in the case of 
[H, M] 0, where H and M are the Hamiltonian and 
the magnetization operators, respectively. Generally, the 
analysis of xv in the interacting system is not easy. 
In the mean-field approximation (Gutzwiller approxima- 
tion), H where we replace Vq with ^fz ■ Vo (z <C 1 is the 
renormalization constant) , both xp and xv are enhanced 
by 1/z. However, this result depends strongly on the 
assumption that the /-electron level Ef is renormalized 
by 1/z to E*±, which is never trivial and should be con- 
firmed in detail. For U > 0, it is also nontrivial 
whether the quasi-particle like spectrum at remains 
coherent or vanishes to become incoherent. 

The aim of this work is to study the electronic struc- 
ture and the magnetic susceptibility of J = 5/2 PAM 
beyond the mean-field theory. Here, we perform the nu- 
merical calculation by using the self-consistent second- 
order perturbation theory (SC-SOPT) with respect to 
(w.r.t.) U. In principle, (SC-)SOPT is valid only for the 
weak-coupling regime. Nonetheless, it is known to be 
very useful for discussing strong coupling properties of, 
e.g., the impurity Anderson model jl(| and the Hubbard 
model. [ p| In addition, we use the local approximation 
(d = oo approximation) , |f2| which is known to work well 
for nondegenerate PAM. |lf| , p^ | We check the higher cor- 
rection terms w.r.t. 1/d for J = 5/2 PAM, and find that 
the local approximation works better due to the orbital 
degeneracy. |Q 

We find that, for the six- fold degenerate case (no CEF 
case) , both xp and xv are enhanced by the enhancement 
factor at the Fermi energy within SC-SOPT. This result 
coincides with that of our analytical study on the basis of 
the Fermi liquid theory for the strong correlation region, 
given in ref. ||. This numerical work is complementary 
to the analytical work in ref. ||. We also confirm the 
renormalization of Ef towards the Fermi energy, /x. This 
result is consistent with the enhancement of xv- 

The Hamiltonian for J — 5/2 PAM is given as 
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where cj_ (/Jm) ^ s ^ ne creation operator of a conduc- 
tion electron (/-electron), and ek is the k-dependent 
spectrum for conduction electrons. M is the eigen- 
value of J z (M = 5/2,3/2,- • -,-5/2). The c-f mix- 
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in g pot ential V k Ma in (§) is given by V kM a = V Q ■ 
V^E™ -cry/(7/2 - Ma)/76 mM _ a/2 Y^(0 k , <p k ), 
where ^™ 3 (^fe,yfe) is the spherical harmonic function. 
We impose the Zeeman energies for only /-electrons, 
E* M = Ef + gpgM ■ H, where g is Lande's ^-factor 
( g = 6/7 for J = 5/2 ), and we set fig = 1 here- 
after. The electronic structure in the case of U = is 
shown schematically in Fig. |l|. A_ = i?fe = 2 W , and E^ 
is the quasi-particle spectrum, which is the bonding or- 
bit brought about by ek and Ef. There remain fourfold 
nonbonding orbits at Ef. 

The local Green's function g(ui), with which the self- 
energy S(w) is constructed in the local approximation, is 
given at H = by §] 

g(w) = ±G(w) + lGf(u), (2) 
G H = Jf E (V^M - ZV 2 /{u + (J, — e^y 1 , (3) 

k 

G f (Lj) = (u + (i-Ef-Z(u)y 1 . (4) 

and the w-dependent renormalization factor in the Fermi 
liquid theory is 

«(w) = (i-as M (w)/aj)- 1 > (5) 

This w-dependence is neglected in the mean-field treat- 
ment. Note that both g(u>) and E(u>) are diagonal and 
independent of M in the case of no CEF. || Considering 
the Pauli principle, the self-energy by SC-SOPT in the 
local approximation is given by 

E(u;) = ^-nf - 5C/ 2 J! ^-g(e)g{e>)g{u + e - e') (6) 

at T = 0. Figure^ snows the calculated self-energy 
The w-dependence of is complicated, and the energy 
region, where renormalization takes place by a factor of 
l/z(u>) > 1, is limited to within 2Tq from the Fermi en- 
ergy. (We call T the coherent energy in the system.) 
These features are not derived by the mean-field treat- 
ment. 

Figure shows the /-electron density-of-states (DOS), 
p(uS) — — ~Im g(u + iS), obtained by using (§)-(§) and 
ml). We also introduce the DOS for the nonbonding or- 
bits, p' (uj) = — -glm G^(uj + iS). In the numerical cal- 
culation, we imposed weak fc-dependence on Ef for con- 
venience. The noteworthy points are as follows, (i) Et 
is renormalized toward p to be Ef, and its weight is re- 
duced. We find that the relation \E*f — p\ ~ z(0) • \Ef — p\ 
holds well for smaller U where SC-SOPT is adequate, (ii) 
Broad satellites (or shoulders) are formed on both sides 
of p, which are attributed to Mott-Hubbard excitation. 
Large part of p> (uj) is transferred to the satellites, which 
is localized at Ef in the case of U = 0. This feature is not 
derived by the mean-field theory. The renormalization of 
Ef observed in our calculation means that Et lies inside 
the coherent region, and \v is expected to be enhanced. 



Next, we calculate the uniform magnetic susceptibility 
for the following two typical cases. 

(i) Case 1 (there is no CEF and three KDs are degen- 
erate): In the previous work based on the Fermi liquid 
theory, M we found the following result in the strong 
coupling region: 

x/x° — 1/ Z ff(0) for the metallic case, (7) 
X/x — l/zff(Al) for the insulating case, (8) 

+ (9) 

whose deviations are of order ~ O(z ), and 1/zb{w) = 
1 + jL^L£m(o;) is the magnetic enhancement factor in 
the Fermi liquid theory. In general, l/z(0) < 1/zh{0) 
through the Landau parameter Fft < 0. j8j Thus, 
Xinsuiator = Xv is enhanced. In the metallic state, total 
susceptibility x is given by the £;-limit of the dynami- 
cal susceptibility XfcM (i.e., x = hm fe ^oXk(0)) and xv 
should be defined by the w-limit of Xfe(w). In the insu- 
lating state, the two limits coincide. 

Using the above relations, we calculate x in a manner 
consistent with the self-energy obtained by SC-SOPT, so 
as to satisfy the conservation law (i.e., the Ward identity) 
proposed in ref. Pq j. Furthermore, we compare it with 
l/z H (0) (or derived using SC-SOPT. Here, 

^7^-5a/( w )|h=o is given by (5.12) in ref. §], and is com- 
plicated because M has interband components. (Note 
that [H,M] 7^ 0.) The results for (a) the metallic case 
and (b) the insulating case at T = are shown in Fig. 
||. We can see that relations (Q) and (||) hold with high 
accuracy for small U, where the result obtained using 
SC-SOPT is reliable. We consider the deviation from (J7]) 
and (||) for large U in Fig. to be artificial because it 
is of the order, at least, 0(U 3 ). Taking account of the 
analitical result in the strong coupling region, || we con- 
clude that the relation (Q) and (§h hold irrespective of the 
strength of the correlation. Moreover, Fig. ||(a) clearly 
shows the enhancement of xv i n the metallic state. 

(ii) Case 2 (where only the lowest KD is active due 
to strong CEF): The system reduces to the familiar non- 
degenerate PAM. In this model, xv is finite unless the 
(/-values for /-electrons (<?/) and for c-electrons (g c ) are 
the same. M Here, we set g c = 0. The result of the 
mean-field approximation is 

XP = 9 2 f ■ p(0)/z, (10) 
Xv = g) ■ Pc ■ cot- x (ViV/(/x - e*,)), (11) 

where p(u) (p c ) is the /(c)-electron DOS. (xp =0 > Xv = ° 
because p(0) = V 2 /(Ef — p) 2 ■ p c .) Thus, the enhance- 
ment of xv is of the order O(l) even in the mean-field 
approximation. This behavior is recognized more accu- 
rately in quantum Monte Carlo simulation within the 
d = oo approximation. Jl6|,[l^] Here, we calculate x an d 
1/zjj for the insulating state using SC-SOPT, and show 
the result in Fig. M. We can see that x deviates from 
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1/zh in the £/ 2 -term. Especially, the sign of curvature of 
X w.r.t. U and that of 1/zh are opposite. 

Here, we compare several properties between case 1 
and case 2. The sign of the curvature of xv w.r.t. U 
are opposite in the two case. Thus, interestingly, orbital 
degeneracy is indispensable for the large enhancement of 
Xv- Next, we discuss the coherent energy To for both 
cases (see Fig.||). We can see that T in case 1 (sixfold 
degenerate) is much larger than that in case 2 (twofold 
degenerate) if 1/z is the same in both cases. Roughly 
speaking, within SC-SOPT, we can see that T /z ~ \Ef — 
n\ for case 1 and Tq/z ~ |A_— fi\ for case 2 in the metallic 
system, where Ef and A_ are not renormalized values. 
This result indicates that the paramagnetic state is more 
stable in case 1. Finally, we discuss the Wilson ratio 
R for both cases, where R = (x/x U= °) ' z m the local 
approximation. In Fig. ^, we plot the ^-dependence of 
R derived using SC-SOPT for both case 1 and case 2. 
In case 1 , R increases much more slowly w.r.t. U . As a 
result, the orbital degeneracy stabilizes the paramagnetic 
state against (ferro-)magnetic instability. This tendency 
is consistent with the discussion on T above, and it is 
also found within the Gutzwiller approximation. 

In this work, we have studied the /-electron DOS and 
X in J = 5/2 PAM by applying the SC-SOPT, in order 
to study many-body effects beyond the mean-field ap- 
proximation. We showed that the Coulomb interaction 
drastically changes the whole electronic structure. This 
change causes the enhancement of xv, which is the in- 
terband contribution away from the Fermi level. In Fig. 
[|, we found that the relations (^) and (||) hold with high 
accuracy, which means both xp and xv are enhanced 
by the magnetic enhancement factor, 1/zh- This result 
coincides with our 1996 analytical study based on the 
Fermi liquid theory. Interestingly, the enhancement of 
Xv is closely related to the orbital degeneracy because 
it is not observed in nondegenerate PAM (see Fig. ||). 
Moreover, comparing the Wilson ratio between J = 5/2 
PAM and nondegenerate PAM, we can see within the 
scope of SC-SOPT that orbital degeneracy stabilizes the 
nonmagnetic state. 
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FIG. 1. Schematic electronic structure of J = 5/2 PAM 
for U = in the metallic state (/i < A_). Fourfold degen- 
eracy remains at Ef. Band insulating state is realized for 
A- < /i < E f . 
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FIG. 2. The self-energy of the metallic J = 5/2 PAM by 
SC-SOPT for U = 0.3, where 1/2(0) = 1.53. For |w| < 0.02, 
the relation ImE A H = U 2 ■ (5tt/2)p(0) 3 ■ oj 2 w 2.4 • uj 2 is 
satisfied. Other parameters are — — 1 + 2- (k/7r)' , Vb = 0.4, 
Ef = —0.25 and /j, = —0.38. In each stage of self-consistent 
calculation, we choose Ef so that fi is unchanged. Note that 
E) w 0.1 < To (see Fig. 3). 
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FIG. 4. In case 1 (sixfold degenerate case): x> V^h and 
1/z obtained using SC-SOPT for (a) the metallic case and for 
(b) the insulating case. For (a), we used the same parameters 
as those in Fig. 2. For (b), parameters are tk — — 1 + 2- (k/w) 3 , 
Vq = 0.4 and E f = 0. We also imposed fi = (A* + E})/2 
which should be satisfied in the case of a pure insulator at 
T = 0. 



FIG. 3. The /-electron DOS of the metallic J = 5/2 PAM 
obtained by SC-SOPT for (a) U = and (b) U = 0.3, where 
l/«(0) = 1.53. Other parameters are the same as in Fig. 2. We 
can see that p^(w) for U = 0.3 is spread into a wider energy 
range. p(0) is unchanged by U in the local approximation. 
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FIG. 5. In case 2 (twofold degenerate case): Xi 1/zh (0) and 
1/2(0) obtained by SC-SOPT for the insulating case. Here 
we used a symmetric PAM (e k = -1 + 2 • (k/71-) 3 , V = 0.4, 
Ef = n = 0). 
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FIG. 6. [/-dependence of the Wilson ratio obtained using 
SC-SOPT for case 1 (i?6-foid) an d for case 2 (i?2-foid)- In both 
cases, we used the same parameters (and p(0)) as those in Fig 
2. 
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